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A few derivation methods of the distribution in the title are discussed from some different 
viewpoints. It is found that this type of fading may be caused by the interferences of two 
correlated waves. A transitional aspect of the distribution due to the change of a parameter 
is illustrated in detail. Some applications of the distribution to signal statistics are dis- 
cussed. The functional similarity is found between the characteristic function and the error 
probability of signals. Using this basic relation some error probabilities are estimated for 
various cases. 

1. Introduction 

Since about 25 years ago, the author lias engaged in an extensive investigation of the 
characteristics of radio fading phenomena. At an earlier stage of our series of observations, a 
kind of fading much severer than the Rayleigh type was sometimes experienced in particular 
circuits, especially in the Taipei-Tokyo circuit operating at the frequency of about 10 Mc/s. 

The appearance of such a type of fading was quite unexpected. Therefore our considerable 
interest was directed to finding the causes of this type of fading. Unfortunately, however, any 
convincing explanation about the mechanisms of the fading seemed to be beyond our knowledge 
at that time. 

In due course, we fortunately had an opportunity to make careful observations on several 
communication circuits in order to find possible mechanisms of lading, using much improved 
equipments and advanced techniques of observation. Through this well-prepared experiment, 
the existence of a type of fading stated above was clearly confirmed. 

Meanwhile, our theoretical study of fading was proceeding a little ahead of our experi- 
mental research. In the course of our statistical approach, three basic distributions w T ere found 
as particular solutions of the so-called problem of random interference. One is the distribution 
referred to in the title, tentatively named the "^-distribution" [Nakagami and Sasaki, 1943]. 
Soon after that, it was found that the formula accounted for the type of observed fading which 
had shown much larger fluctuation than the Rayleigh fading. 

It seems to the author, however, that most of the researchers in this field still have a mis- 
conception that the Rayleigh distribution exhibits the largest fluctuations among all types of 
fading. 

In the following, some discussions will be made on the ^-distribution. In particular, we 
shall derive the distribution on the assumption of some propagational modes which may be 
supposed to exist in actual fading. 

2. Derivation of the q-Distiibution 

We shall derive the ^-distribution and its generalized form from gaussian distribution. 
Now; let us start from the two-dimensional gaussian distribution 

^ H 2wJ 2 a-v) exp {-2oh)^r 2 M + '0}- (1) 
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By changing the variables (x, y) to (R cos 6, R sin 0), and performing the integration with 
respect to 6 from to 2ir, we have 



p(R)= *p(R cos0, RsmS) 



b(x,y) 



2>(fi,< 



^0 



V<Ti(T 2 (l — P 2 ) 

Here, if we put 



e-W^lK, '— V7(r-f: 



^ + M^VKH)+(vfe)> 

f (3) 

then (2) takes a standard form of the g-distribution 

It is of great importance to note that either in case of a 1 equal to cr 2 , or p 2 equal to zero, (2) 
retains the same standard form as (4). However, if the above two conditions are satisfied 
simultaneously, the distribution turns out to become the Rayleigh form 

p (R) = — : — e~^+w. (5) 

0"l + 0"2 

Further, if p tends to unity, then approaches to zero, and (4) takes the following form 

1 R2 

p(R)=- 1 ===e <n+*2. (6) 

It is of much interest for us to see that the limiting form (6) is identical with the limiting case 
of the observed m-distribution [Nakagami, 1943], 

r>J?2m-l ry) m m R 2 1 

P(fl)= r(nOCr '~-m(B,m,Q), m>^ (7) 

as the parameter m approaches }{. In this regard, we shall discuss the fading mechanisms 
further in the following section. 

2.1. Derivations of a Generalized and an Additive Form of the qr -Distribution 

A generalized and an additive form of the distribution may be deduced in several ways. 
The simplest cases, however, are shown as follows: 
First, if we let 

1^12^-1 _j2 I |2m-l _p 

P( X ) = Lw \ e a >V(y) = em-nt \ e * > ( 8 ) 

1 a m T(m) ><r^»' (3 m T(m) 

then we may have the standard form 

PW- n 1V n-i /2 ^-i/2i yU--) r (9) 

r(ro)(a0)»Q--±J ^ V/ * a/J 

The detailed properties of this type of distribution are fully discussed in Nakagami and Nishio 
[1954a]. 
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Next, we shall derive a simplest additive form of distribution. Now we start with the 
distribution 

p{x>y)= (*y) 2 exp f _*=*S3+l£\. (10) 

* 2^<t 3 (1+2 P 2 )Vp(1-P 2 ) I 2(i- P 2 y / v ; 

By the transformation of variables from (x, y) to (R cos 6, R sin 0), we have the distribution 

where 

a=2<j(l + p),P=2a(l-p). (12) 

3. Derivation of the q-Distribution From an Actual Mode of Propagation 

We are now in a position to derive the ^-distribution from the two correlated Rayleigh 
or the ra-type of fading. At first we shall begin with the Rayleigh type of fading in the following. 

3.1. Case of Correlated Rayleigh Fading 
As is well known, the joint distribution of the Rayleigh type assumes the form 

Qsrf+Oirj 



4/y>2 



*<*. *-ssfc*j - '""■"-'" '- tera- (13> 

where U l = ri and fi 2 =/ , l, respectively, and k 2 denotes the correlation coefficient between 
r? and v\. 

Then, the resultant intensity distribution of the sum of two vectors f x and r 2 and a scat- 
tered component can be easily expressed by the Hankel form of integral: 



V 



(K) =B r\J (\R)F(\)dk, (14) 



where F(X) is tentatively called the amplitude characteristic function. If we assume that the 
phase difference of f x and r 2 is random, then F(\) becomes 



F(\) = J (\r l )Jo^r 2 )e 4 , (15) 

where 5 indicates the scattered power, and 

Jo(\r 1 )J (\r 2 )= j I MXrO Jo(\r 2 )p(r h r 2 )d ri dr 2 . (16) 

o 

This form of integral may be readily evaluated as 

Jo(Xr,)Jo(Xr 2 ) = e^ ( " 1+ " 2) / (S^. (17) 

Hence, we may arrive at final result 

P (R) =B (\j„(X7?)/„ (S^) <f* «**« dx 

= v^T ■^'Mw-*)/* (18) 
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where 



a=(Qi+0 2 +5)+2VQi0 2 t,l 



(19) 



j8*=(0 1 +Q f +*)- 

In a limiting case of Z: equal to zero, ^?(i?) assumes the Rayleigh form 

2R 



P(R) = 



e J2i+0 2 +5. 



fii+122+5 
Under the conditions ftj^fta^ft, k-^1, and <5-^>0, (18) turns to the half Gaussian form 



p(R) = 



1 -£ 



(20) 



(21) 



It is of great interest to see (20) and (21), especially the latter distribution; because these 
two distributions are both situated at the extreme boundaries of the ^-distribution, as was 
already proved from a different viewpoint. 

In figure 1, a family of curves of the distribution are shown with e = 8/tt as a parameter. 
They are very suggestive of the transitional aspect of the shape of distribution. 

The comparison of the ^-distribution with the m-distribution is shown in figure 2 under 
the extreme conditions respectively. From this figure we may find a good agreement of the 
two formulas even at the extremities. This functional similarity was already discussed in 
detail. 

Some examples of the deepest type of fading observed in a microwave communication 
circuit in Japan are shown in figure 3. The foregoing discussions seem to maintain that the 
observed distribution shown in figure 3 might be caused by correlated interferences. 
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Figure 1. Transitional aspect of the shape of the Figure 2. Comparison of the q-distribution with the 

probability density function p(x). m-distribution having the parameter m close to }i. 
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Figure 3. Comparison of the observed distribution in 
a microwave circuit with the m-distribution having 
the parameter m smaller than 1. 



20 
10 





































4 




Date; Augusf 5, (954 

Frequency; 6 Gc/s 

Distance; 1 50 km 

f-o- oo.oo-oi.oo 

T ' m «; <^ 04.00-05.00 
(J ST) [_ M _ 06 . 00 _ 07 . oo 










) 


If 








I 








fl 
































1 
















I 
















V 






























/$' 
















f 





























'm«Q7 








m = 


05 \ 




^ 








m = 



























-60 -50 -40 -30 -20 -10 <0 20 

Re/ative Intensify in dB 



3.2. Derivation From the Correlated m -Distribution 

We now proceed to a more general case of two correlated m-variables. As is well known, 
the joint probability density function assumes the form [Nakagami and Nishio, 1954b, 1955] 



p m (r u r 2 ) = 



Q2Tl-\-<T\T2 

4(r 1 r 2 ) m e _ (n «(i-(i J ) 



r(m)a 1( r 2 (l-F)(V^^)" 



-/„_,{ »™ \, (22) 



where 

o- im =ri, o- 2 m=rl. Therefore 



where 



p(R)=R\ XJo(X^)^o(Xr 1 )J (X^)^ d\ 

M\r 1 )J Q (\r 2 ) = J J Jo(Xr 1 )J (Xr 2 )p m (r 1 ,r 2 )^r 1 rfr 2 



(23) 



(24) 



For the sake of simplicity, we treat the simpler case of m=2 and 5=0. The calculation 
of (23) is somewhat complicated even under the simplified conditions, a l = o- 2 and m=2, 
respectively. The final result is 






(25) 



This form of distribution resembles the distribution (11). In this case, if k 2 tends to zero, 
(25) takes the form 

p(B)=h{ 771 (R, 1, Q)+ 771 (R, 3, 12) }. (26) 
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4. Some Applications of the Distributions to Signal Statistics 
4.1. Criterion of Signal Improvement 

In the field of signal statistics, it is important to evaluate the performance improvement 
attained in various communication systems. For this purpose, two criteria may be used. 
One is expressed by the cumulative distribution, P(R), and the other by the so-called error 
probability, P e . 

The former is expressed as 

1 (*c + j«> zR 2 J„ 

P{R)~A e *( 2 )f> (27) 

(*™ - zR 2 

4>{z) = e p(R)dR, (28) 

where P(R) and p(R) are the probability and the probability density function of signal in- 
tensity R. This criterion apparently seems to be irrelevant to noise information. The latter 
criterion, P e , however, contains noise information explicitly. Though it is a quantity defined 
with respect to digital communication systems, it can be used in other systems as well. As 
we are here concerned with fading signals, we may better choose the average error probability 
P e as an adequate measure of signal quality. 

In general, P c takes various forms depending on the methods of signal detections. In 
the case of unknown phase information [Turin, 19581, it is simply expressed by 

P=± e -^> (29) 

where E is the signal energy and N the white noise power density, respectively. The average 
error probability, P e in this case, is simply expressed by 






e~*Np(R)dR (30) 



=J[*(2)U' (3D 

where E=\R 2 T ) T— element length in seconds. 

It is worthy of notice that P e assumes the identical form with the characteristic function, 

Through the author's experience, it may be said that (f>(z) usually takes a simpler form than 
p(R). This fact will be of great help when the treatment oip(R) is complicated. 

On the other hand, it would be more convenient, in many cases, to use the Hankel type of 
integral rather than the Laplace transformation. However, as is mentioned above, <$>(z) is 
usually more suitable to estimate P e . Therefore the following well-known transformation 
from F{\) to <j>{z) 



* w =sf 



e"zF(\)d\ (32) 



will be useful. 
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4.2. Error Probability 

For the estimation of average error probability, as given by (31), it is sufficient to find the 
characteristic function, <t>(z). On this basis, we shall proceed to calculate the error probability 
of combined signals in diversity reception. 

a. Case of Squared Addition 

Let us begin with the case of the sum of squared signals: 

R 2 =Rl+Rl+. . .+R 2 n . (33) 

The characteristic function, <j)(z), is expressed by 

*(*)= J°e-«** +I * + - ■ - +Rl) p (fi u R 2 , . . ., R n ) dR t dR2 ■ • . dB„ (34) 

where p(Ri, R 2 , . . •, R n ) is the joint probability density function of /(-diversity signals. 

For simplicity, we assume that the signals are independent of each other, then (34) takes 
the form 

0(s)=n0,(z). (35) 

i=\ 

Accordingly, P e is simply expressed as 

P e = 2^U i P e) (36) 

i=l 

where t P e shows the error probability of the it\\ branch of signal. 

Further, we shall consider the case of two correlated ra-fading signals; then 



00 

4>{z)= J*JV ,t, * fi *fiw(fii, B 2 )dR l dR2, (37) 

o' 



where p m {R\, 7i%) stands for (22). After some calculations, P e is obtained as 



Pe 2 [_{ l + 3(<r 1 + <r 2 )+ 2 2 (l-F)<r 2 «T 2 } m l=-£ 



b. Case of Linear Addition 



(38) 



We shall treat the sum of two independent signals which follow t lie /^-distribution. In 
this case, <j>(z) takes the following form 



oo 

4>(z)=-- J*JV*«i+«2> 2 777 (R u m> q) ffl (R 2! m, Q) dBtdB a . 



(39) 



Applying the well-known transformation 

Ri=ue- x ,R 2 2 =ue x , (40) 



(39) is brought to 



_4Z^J* V- 1 ^ | 2 (^ +2 j M y e -^du, (41) 



r 2 (m)S2 2m 

where K (x) is the modified Bessel function of second kind. 
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In table lare shown some of the calculated results. 

Table 1 



m 


</>(*) 


1 
2 


2 1 f 2122 
, cos -1 1 Oo 
7r Vl+4i22 1 + 2122 


1 


1 [~1 " 2 co^ Qz 1 




1 + 2122 L Vl+2122 l+«d 




2 


16(lil22)3[ 16 + 16 ^ + 1 ^ 2 

3122(12 + 12122+ 512 2 2 2 ) 

; COS" 

2V1 + 122 


1 ® Z "1 

2+Q2J 



c. Case of Vectorial Addition 
Next, we shall proceed to a more general case of two vectorial sum 

R 2 =Rt+Rl+2R l E 2 cos e, 



(42) 



where we assume that R x and R 2 follow the correlated m-distribution and varies at random. 
Then 



LIT oo 

«*-hStt 



-z(R\+Rl-\-2RiR2eosd) t^n7^» 

e Pm(Rii R2)dBidR 2 d9 







which yields, after some calculations. 

4 



*(») = 



r(m)cr 2 (l-P)( 



>*) m Jo 



^ w i^o{2(a+s)^}/ w _i(2(2to)7 (2^)^, 



(43) 



(44) 



where 



a 1 =cr 2 =o- J a- 



1 



(!-#> 



For larger values of m, the calculations become more complex. The results are shown in 
table 2, where calculations are limited to smaller values of m. 





Table 2 


m 


0(g) 


1 


1 




V{ 2<r2(l+fc) + 1 H 2(T2(1- A;) + 1 ) 


2 


2(T2J(r(l-/b 2 )2+l( + l 


[{ 2az(l-k) + 1 J { 2<r*(l + fc) + 1 |j 3/2 




1 + 4(T2{ 1+CT2(1-/C 2 ) | + 6<x 2 > 2 { 1+<T2(1-A; 2 ) j 2 


3 


[l + 4(T2 + 4<7 2 2 2 (l-/r 2 )] 5/2 
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